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Abstract. We consider the continuous version of the Vicsek model with noise, proposed 
as a model for collective behavior of individuals with a fixed speed. We rigorously derive 
the kinetic mean-field partial differential equation satisfied when the number A*' of particles 
tends to infinity, quantifying the convergence of the law of one particle to the solution of the 
PDE. For this we adapt a classical coupling argument to the present case in which both the 
particle system and the PDE are defined on a surface rather than on the whole space R'*. As 
part of the study we give existence and uniqueness results for both the particle system and 
the PDE. 

Introduction 

The stochastic Vicsek model [13j arises in the study of collective motion of animals and it 
is receiving lots of attention due to the appearance of a phase transition [21 [9] . A continuum 
version and variants of this model have been proposed in the recent works OH]. Our objective 
is to rigorously derive some continuum partial differential equations analysed in [5] from the 
stochastic Vicsek particle model. This was carried out for a family of collective behaviour 
models in [1] following the method of [12]. The present models do not fall into this analysis 
due to the evolution being defined on a surface as we explain next. In the models considered 
here, individuals are assumed to move with a fixed cruising speed trying to average their 
orientations with other individuals in the swarm in the presence of noise. This orientation 
mechanism is modelled by locally averaging in space their relative velocity to other individuals. 
More precisely, we are interested in the behaviour of N interacting M^'^-valued processes 
{XI, V^)t>o with 1 < i < with constant speed |VJ*|, say unity. We define them as solutions 
to the coupled Stratonovich stochastic differential equations 

(dxi = v^dt, 

dV; = ^/2 o dBl - P(yi) I ^ > : K{Xl-Xim -y^)\ dt. 



V2P{Vl)odBl-P{Vi) \-^Yl^K{Xl 



Here P{v) is the projection operator on the tangent space at v/\v\ to the unit sphere S in 
W^, i.e., 

V ®v 



P{v) 



This stochastic system is considered with independent and commonly distributed initial data 
{Xi, V^) eW^x S"'-^ with l<i< N. The {Bl)t>Q denote N independent standard Brownian 
motions in W^. The projection operator ensures that V^* keeps constant norm, equal to 1. The 
second term in the evolution of VI models the tendency of the particle i to have the same 
orientation as the other particles, in a way weighted by the interaction kernel K, as in the 
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model proposed by F. Cucker and S. Smale [3]. Let us observe that P{V^)V^ = 0, so we can 
drop the corresponding term when writing ([T]) to recover the usual formulations as in j5]. 

We will work with stochastic processes defined on R^rf instead of M*^ x S'^-^ We will check 
later on that solutions of ([1]) with initial data in R'^ x S"^"^ remain there for all times. We have 
written ([T|) in the Stratonovich sense, since the term involving noise corresponds to Brownian 
motion on the sphere S*^"^ as in [TO', Section 1.4] and [HI Section V.31]. 

By symmetry of the initial configuration and of the evolution, all particles have the same 
distribution. Even though they are initially independent, correlation builds up in time due to 
the interaction term. Nevertheless, this interaction term is of order and thus, it seems 
reasonable that two of these interacting particles (or a fixed number k of them) become less 
and less correlated as N gets large (propagation of chaos). 

Following [12] we shall show that the N interacting processes {X'l:,V^)t>o respectively be- 
have as — 7- oo like the auxiliary processes {Xl,Vl)t>o, solutions to 

< dy; =y2Pivi) o dBi - p{r,)iH^* mxi, v\) dt, (2) 
^(x^,Fi,) = (x^,yj), /i = iaw(x;,r,) 

in the Stratonovich sense. Here the Brownian motions {Bl)t>Q are those governing the evo- 
lution of the {Xl^ Vl)t>o and 

{H*f){x,v)= [ K{x-x'){v-v')f{x',v')dx'dv', x,v£M'^. 

Note that ([2]) consists of N equations which can be solved independently of each other. Each 

of them involves the condition that ft is the distribution of (X^, Vl), thus making it nonlinear. 

The processes {Xl, Vl)t>o with i > 1 are independent since the initial conditions and driving 
Brownian motions are independent. 

We will show that these processes defined on M^'^ are identically distributed, take values 
in X §'^~^ if initially so, and their common law ft at time t, as a measure on M'^ x S'^"^, 
evolves according to 

dtft + io-V^Jt = A^ft + V^-{ft{I-io(^io){H*ft)), t >0,xGM^l^gS'^-\ (3) 
Now the convolution H * f is over M'' x 

{H*f){x,uj)= [ K{x-x'){uj-oj')f{x',J)dx'dJ, x G M^, w G S'^"^ 

Moreover, V^; stands for the gradient with respect to the position variable x G M"' whereas V^j, 
Vo;- and respectively stand for the gradient, divergence and Laplace-Beltrami operators 
with respect to the velocity variable a; G S'^~^. 

This equation is proposed in [4J as a continuous version of the original Vicsek model |13j . 
and one of our purposes is to make this derivation rigorous. The asymptotic behavior and 
the appearance of a phase transition in the space-homogeneous version of ([3]) (i.e., without 
the space variable) has been recently studied in [7]. It is also known as the Doi-Onsager 
equation, introduced by Doi in [6] as a model for the non-equilibrium Statistical Mechanics of 
a suspension of polymers in which their spatial orientation (given by the parameter oj G S"'"^) 
is taken into account. 

The main result of this paper can be summarized as: 



MEAN-FIELD LIMIT FOR THE STOCHASTIC VICSEK MODEL 



3 



Theorem 1. Let /g be a probability measure on R"' x S^^^ with finite second moment in 
X S M*^ and let [Xq, Vq) for 1 < i < N be N independent variables with law /q. Let also K be 
a Lipschitz and bounded map on W^. Then, 

i) There exists a pathwise unique global solution to the SDE system ([T]) with initial data 
(Xq,1q*) for l<i< N; moreover, the solution is such that all have norm 1. 

ii) There exists a pathwise unique global solution to the nonlinear SDE ([2]) with initial 
datum (X'q,Vq); moreover, the solution is such that V\ has norm 1. 

iii) There exists a unique global weak solution to the nonlinear PDF ([3|) with initial datum 
/o- Moreover, it is the law of the solution to ([2]). 

Solutions to general SDE's can be built in submanifolds of by means of the Brownian 
motion of the ambient space as in [11^ Theorem V.34.86] for instance; then one can interpret 
the generator in terms of the corresponding Laplace-Beltrami operator. For example, the 
Brownian motion on a submanifold S of is the solution to the SDE 

dWt = P^{Wt) o dBt 

on and with Pj]{w) being the orthogonal projection of onto the tangent space at w to 
E. Here, we give the full construction and derivation of the evolution of the law as it can 
be done explicitly in the case of the sphere S"^"^. Let us also emphasize that we have only 
partial diffusion since it is a kinetic model. 

We observe that existence of L^ and classical solutions for the space-homogeneous version 
of ^ has also been considered in [7]. 

As a direct consequence of the classical Sznitman's theory, we get the following mean-field 
limit result: 

Theorem 2. With the assumptions of Theorem{l\ and for the respective solutions {XI, V^)t>o 
and {Xl,Vl)t>o of ([T]) and for allT > there exists a constant C such that 



E 



q-X\^ + \vi-v\\^ 



c 

< — 

- N 



for alio <t <T, N >l andl<i<N. 

This estimate classically ensures quantitative estimates on (see ||12i, Qj for details) 

i) the convergence in of the law at time t of any (by symmetry) of the processes 
{XI, Vl) towards /t, 

ii) the propagation of chaos for the particle system through the convergence of the law 
at time t of any k particles towards the tensor product fP^ (for k fixed or k = o{N)), 

iii) the convergence of the empirical measure at time t of the particle system towards ft- 

Of course, the same techniques lead to a corresponding mean-field limit result for the space- 
homogeneous particle system instead of ([T]) , obtaining the corresponding space- homogeneous 
PDE. 
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Proofs 

Using the standard Ito-Stratonovich calculus, see [U p. 99] for instance, equations ([I]) 
and ([2|) are respectively equivalent to the Ito stochastic differential equations 



'dxi = vidt, 



N 



dy: = V2P(vi)dBi-p{vi) I -Y^K{xi-xi){vi-v/) I dt-{d-l) 



dt. 



and 



dXt = Vtdt, 



dt. 



(4) 



(5) 



dV\ = V2P{V\)dBl - P{V\){H * /i)(X;,r,) dt-{d- 1) 

{%Xo) = {XlVi). /i = iaw(x;,rj 

which we now consider. 

We start with the proof of Theorem [H We use a regularization of the diffusion and drift 
coefficients. We let cti be a d x d matrix valued map on with bounded derivatives of all 
orders such that ai{v) = P{v) for all v with \v\ > 1/2, and cj2 and 173 be maps on R'^, again 
with bounded derivatives of all orders, such that cj2{v) = v/\v\'^ if \v\ > 1/2 and cr3{v) = v if 
\v\ < 2. 

Existence and uniqueness for the particle system ([4]). Given such ai,a2, the system 
of equations 

'dxi = v;dt, 

dVi = V2a,{Vi) dBi - a,{Vi) Jl^ K{Xl-Xl){Vl - F/) jdt - {d - l)a2{Vl)dt 

(6) 

starting from (Xq, V[{) G M*^ x ^ for 1 < i < has locally Lipschitz coefficients. Moreover, 
by the Ito formula and as long as |V^'| > 1/2, 



d\V 



i\2 



■ P{V')dB' -2V'- P{V') ^^Y1k{X' - X^){V' - W)j dt 



d 



-2id-l)dt + 2Y,5,,d{Bl-Y,^Bi,,Bi-Y,jL^Bl 

q=l 



kl=l 



-2 (d - l)dt + 2 ^ 



k=l 



\ p=l 

p=i 



w 



i\4 



dt = 0. 



Here we dropped the time dependence, wrote y = (yi, . . .yd) £ and used the fact that 
■ P{V^)y = for all vectors y G W^. Hence \V^ \ = 1 up to explosion time. Since moreover 
dXl = VI dt, this ensures that the explosion time is infinite, hence global existence and 
pathwise uniqueness for 
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Now the solution to Q for given cJi, o"2 is a solution to ^ since all velocities have norm 1, 
which provides global existence of solutions to (jll). If now we consider two solutions to ([4]) 
for the same initial data and Brownian motions, then they have velocities equal to 1, so that 
are solutions to ^ for any ai,a2, for which pathwise uniqueness holds: hence they are equal. 
This proves the first part in Theorem [TJ 

Existence and uniqueness for the artificial processes ([5]). Let <Ti,a2,(73 be any maps 
as above and let 

Ha-AfKx) = / K{x-y)a3{v-w)f{y,w)dydw. 
Jm^d 

Then, given a distribution /o on M'^x S"^ ^ with finite second moment in x G M.'^ and {Xq, Vq) 
with law /o, the nonlinear equation 

'dXt = Vtdt, 

dVt = V2ai{Vt)dBt - ai(yt){H^, * ft)(Xt,Vt)dt - {d - l)a2{V t)dt, (7) 

has bounded and Lipschitz coefficients on M^"', so admits a pathwise unique global solution 
according to [121 Theorem 1.1]. Moreover, as long as \Vt\ > 1/2, then we can repeat the 
argument above to prove that (i|T^tp = 0, so that \Vt\ = 1 for all time. In particular 
the obtained solution {Xt.,Vt)t>o is a global solution to the genuine nonlinear equation ([5]). 
Pathwise uniqueness of solutions to ([5]) can be obtained as for ^ . 

Existence and uniqueness for the PDE ([3]). Let /o be a distribution on M*^ X S'^-i with 
finite second moment in x S M.'^, {Xo,Vo) with law /o, and let {Xt,Vt)t>o be the solution 
to ([5]) with initial datum {Xq,Vo). Its law ft, as a measure on M?'^, satisfies 

^ / (pdft= [ (v-Vx^ + UesSyip : {I-v0v) + V^ip-{I-v(S)v){H*ft)-{d-l)v-Vy(p)dft 

dt JK2d J^2d 

for all smooth ip on R^*^ by the Ito formula; here and A^, are respectively the gradient and 
Laplace operators with respect to u G W^, and Hess^c/? : M is the term by term product of 
the Hessian with respect to v matrix of (p with a matrix M. 

We have observed that \Vt\ = 1 a.s., so ft is concentrated on M"' X S'^ ^. We now define 
the restriction Ft of ft on M"^ X S''-^ by 



/ ^dFt= [ 



fdft 



for all continuous maps $ on M'^ X S"^ ^, where ip is any continuous and bounded map on R^"^ 
equal to <I> on M"^ x S'^^^. Let now ^ he a. map on R'^ x S'^^^ and ip he a map on M^"^ 
such that <p{x,v) = ^{x,v/\v\) for all 1/2 < \v\ < 2. Then ip is 0-homogeneous in v in the 
annulus 1/2 < \v\ < 2, so that v ■ V„(/3 = for all {x,v) in the support of ft- In particular 

^[ <^dFt = ^ I ipdft= [ (vV^ip + A^ip + V^ip-{I-v(^v){H*ft))dft. 

Then the maps v ■ and v ■ ^x'P are equal on M*^ x S'^^-'^ since <I> and ip have the same 
x-dependence. Moreover, V^j^ = Vytp and A^<1> = A^(p> for {x,ui) G M*^ x S*^"^. This last 
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point can be checked by direct computations. Hence 
d 



dt 



[ ^dFt= [ {u}-V^^ + A^<i> + V^^-{I-u}0uj){H*Ft))dFt. 



This ensures that Ft is a weak solution to ([3]). 

We now turn to uniqueness of solutions to For that purpose we let and be 
two solutions with the same initial datum /o, and at each time t we view them as measures 
on M^'^ concentrated on the surface R'^ x We let {xl,vl)t>o and {xf,v'^f.)t>o be the 

solutions to d?]) with drift given by Hfj^ * fl and i/0-3 * ft respectively, and common initial 
datum {Xq^Vq) with law /q. Then their respective laws g} and gl, as measures on R^'^, are 
solutions to the linear PDE 

dtg\ + V ■ V^gl = {{aiaDki g[) + V, • [gl [a^ {H^, * ft) + {d - 1)02)] . 

Since fl is also a measure solution to this linear PDE on M^"^ with bounded and regular 
coefficients, for which uniqueness classically holds, it follows that g] = fl {i = 1,2). Con- 
sequently, the {X\^V^^t>Q are solutions to the nonlinear SDE d?]), for which we have al- 

~ — 1 — 1 — 2 — 2 

ready proved uniqueness. Hence {XfiVt)t>o and {Xt,Vt)t>o are equal, and in particular 

fK=9l) = {9^=)fl 

Proof of Theorem [21 Since \Vt \ = \V\\ = 1 for all i and i, the processes {Xt,Vt)t>o and 
QCfi V\)t>o are solutions of the corresponding equations with bounded and Lipschitz diffusion 
and drifts coefficients as in d?]). Hence we may apply the estimates in [12'j Theorem 1.4] to 
obtain Theorem [2j 
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